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Abstract 

A possible formation of bound state is predicted in the Standard Model extension with with 
additional scalar singlet. A suitable method of solution of the Bethe-Salpeter equation for a 
Higgsonium bound state is proposed. With the help of an integral representation the results are 
directly obtained in Minkowski space. The appearance of Higgsonia is shown to be quite natural 
for considered extension of Standard model and the calculations are presented for the bound state 
composed from two scalars which result from the large mixture of electroweak doublet and singlet 
eigenstate. 
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I. INTRODUCTION 

The Standard Model (SM) of particle physics involves very minimalistic idea of elec- 
troweak symmetry breaking scenario. Consequently SM results with the single real scalar 
field representing the last experimentally unobserved particle of SM - the Higgs boson. 
However, being inspired by the family pattern of SM fermionic sector, it is quite natural to 
consider an extensions of SM with more rich structure of the scalar sector. By doubling the 
doublets or/and adding the gauge singlets, such next-minimal extensions of SM model have 
been considered and studied from various perspectives. Clear motivation for such extensions 
is to reduce some SM shortcomings, like the better agreement with precision electroweak fit, 
the theoretical problem of mass hierarchy or the dark matter problem. Adding the singlet 
scalar is worthstanding in recently constrained SUSY models as well as 

Having more interacting scalar bosons, one can expect qualitative changes in the scalar 
boson sector. In some circumstances the binding forces between scalars can appear strong 
enough to produce bound states. What is the spectrum of appearing bound states and how 
they exhibit their existence in the collider experiments are important questions. To address 
this issue is in general a nonperturbative task and one has to make a reliable approximation. 
The later should be principally improvable to check the consistency of the assumptions made. 
To find a spectrum we assume that bound states are sufficiently stable to be decoupled from 
the continuum so they can be identified by the solution of homogeneous Bethe-Salpeter 
equation (BSE). 

To avoid extraordinary heavy states within strong couplings a specific Higgs selfinterac- 
tion and extension of Higgs sector is required 0]. As we know the light Higgsonia could 
appear when more interacting scalars are present in the nature. So observing a Higgsonium 
resonance at LHC could be another evidence of the physics beyond the SM, whose details 
remains to be studied. As a simplest model we choose for our actual calculation is the 
extension of SM that involves the addition of a real scalar singlet S to the SM Lagrangian. 
The phenomenological implications for singlet extension SM (xSM) has been studied from 
the collider and cosmological perspectives [3|. The later typically require small mixing with 
the SM Higgs and from the perspective of bound states it reduces to the SM. In such cir- 
cumstances it was shown in j^] that super heavy Higgs m# ~ ITeV would be needed to 
form a bound state. In our model we will consider large mixing, which gives us two scalars 
Hi and H 2 both having masses at few hundred GeV. Such model is a slight generalization 
(without Z2 symmetry) of universal doublet-singlet Higgs couplings model |4j. Assuming 
that trilinear couplings have a little effect on electroweak precision observables S and T then 
the lighter scalar is experimentally restricted by the SM Higgs alone, i.e. its allowed mass 
is between 114GeV and 145GeV, while preferred value of heavier H 2 should lie in the range 
~ 200 — 250GeV 0]. In such circumstances, the most promising candidate for Higgsonium 
consist of two H 2 with a mass being not far form the sum of their on-shell masses. The main 
motivation of this paper is to find the explicit solution and figure out the conditions under 
which the Higgsonium is formed. 

II. XSM WITH LARGE MIXING OF SCALARS 

The Lagrange density for the xSM model is 

C = (D^D^H + ^Sd^S - V(H, S), (2.1) 
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where H denotes the complex Higgs doublet and S the real scalar. The term linear in S is 
chosen to vanish after the spontaneous breaking. The potential is given by 

2 r 

V(H,S) = X(H^H - l —f + 6 -±H ] H S (2.2) 

+ 5 1 H \ H S 2 + 5 lV 2 S + ^S 2 + ^S 3 + 
2 2 3 4 

In unitary gauge the charged component of the Higgs doublet H becomes the longitudinal 
component of the charged iy-bosons and the imaginary part of the neutral component 
becomes the longitudinal component of the Z-boson. In unitary gauge the Higgs field doublet 
reads 

Within our notation the mass terms in the scalar potential become 



(2.3) 



Knass = \ (/ h h 2 + + fighS) , (2.4) 

where 

Hi = 2Xv 2 

Z-4 = k 2 + S 2 v 2 /2 

tils = 5 lV . (2.5) 

The mass eigenstate fields Hi 2 are linear combinations of the Higgs scalar field h and the 
singlet scalar field S. Explicitly, the inverse transformation reads 

h = c Hi — s H 2 ; 

5 = s Hi + cH 2 , (2.6) 
where c = cos8, s = sinO and the mixing is determined as 

tanfl = j -, x = ^^ T . (2.7) 

for x > 1. For a heavier singlet (x < 0) we have the mixing angle 



x 
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tan 9 = ^-j . (2.8) 



x. 



Note, here we differ from the convention of @, Q by the sign of , 

The small Si mixing scenario with the light singlet like scalar has been considered in Q. 
For such a case we do not expect bound state since the triplet interaction between scalars 
is very weak. The terms in the scalar potential that break the discrete S — ► — S symmetry 
are proportional to the couplings Si and k 3 and we do not consider these terms as very 
small ones, instead we assume they are large enough to make a bound state providing also 
sufficiently strong communication with the rest of the SM. 

In this paper we will consider relatively large mixing angle 9 with Higgs masses close 
to few hundred GeV which could be promising for experimental observation thorough the 
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LHC era. For such case the constraints from electroweak precision observables and their 
implications for the LHC Higgs phenomenology has been already analyzed in [3| . 
The mass eigenstates are in any case 



Ml 



Ac 2 + /4 s2 + VhS cs 



2 2, 2 2 
fl h 8 + fi s C 



L 2 — h*h° "T" PS"- ~~ H<hS^° ■ (2-9) 

The rest of the scalar potential contains four new parameters which are added to the SM. 
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(2.10) 



Based on the nonrelativistic consideration, the trilinear interaction should be sufficiently 
enhanced against the quartic one, otherwise the bound state cannot be formed. 

The purely cubic interaction between mass eigenstates can be written in the following 

way 

Vcu* = gmHf + <?ii2#i#2 + <7i22#i# 2 2 + 9222HI (2.11) 
with the new couplings related to the old ones as 
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(2.12) 



The quartic couplings between the physical states are related to the original Lagrangian 
parameters in the following way 
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III. HIGGSONIUM IN QUANTUM FIELD THEORY 

Higgsonium, being a two body bound states of two scalars is described by two body BSE: 



VG®T 



(3.1) 



where we use shorthand notation J k = i J anc [ qM j s the two particle propagator of 
the constituent Higgsies, just for now let us consider Higgsonium made of two Hxs. In 
momentum space it can be conventionally written as 



G [2] (k,P) = D(k + P/2, Ml)D{-k + P/2, M 2 ) ; 



(3.2) 
(3.3) 
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Let assume that the attractive interaction between heavy Higgsies H± is strong enough 
to form a bound state. Within the xSM the irreducible BSE kernel in lowest order reads 



V = GXuu 



E 

x=s,t,u 



" 4gm 

X - Ml 



gggm 
Ml 



x 



(3.4) 



where the first term represents a pure constant interaction and the s, t, u are the usual 
Mandelstam variables. 

In the next it is advantageous to explicitly divide the solution which is independent on 
the relative momentum of the constituents. Following the notation in [5| the original BSE 
can be rewritten in the following form 



T p (p,P) = T I (P) [ V p (k,p,P)G®(k,P) + I V p {k, V ,P)G [2 \k,P)T p {k,P) 
j k j k 



where 



Vi = v c + v a , v p = v t + v u 



(3.5) 



(3.6) 



(to avoid confusion with the particle content, here we use the letter / instead of index 1 
and p instead of 2 originally used in the paper (EJ) The first term collects all the constant 
term, i.e the ones that do not depend on the relative momentum. In our tree level kernel 
approximation it reads 



Vj — 6A2222 + 



+ 



36^22 



P 2 -Ml P 2 -Ml 



So the full solution of BSE is given by the sum 

r(p,p) = T I (p) + T p { P ,p), 

where the equation for the function Tj(P) is purely algebraic 



IMP) 



VjJ k T p (k,P)G^(k,P) 



(3.7) 



(3.8) 



(3.9) 



The BSE represents the singular equation which can be solved by some known method. 
One known possibility is to perform a Wick rotation for a relative momenta of constituents 
while keeping the total square of four momenta P 2 timelike. 

The other well known possibility is the Minkowski solution performed within the utiliza- 
tion of the known integral representation for the kernels and amplitudes that appear in the 
BSE. Then the momentum space BSE can be rewritten to the equivalent integral equation 
for the weight function. This, historically named Perturbation Theory Integral Representa- 
tion [s| has been successfully utilized to solve relativistic bound state problems in various 

SB 



models and approximations 



12 . 13| . We straightforwardly apply the method here. 



The bound state vertex function can be expressed as 



1 CXJ 

T(P,p) = J dr] J da 



[F(a, V] P,p)Y 



(3.10) 
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where n is an arbitrary integer and all the singularities are trapped by the zeros of the 
denominator in (13.101) which reads 



p2 

F(a,-n; P,p) = a - (p 2 + P.pz + — ) - ie (3.11) 

Here, having the quartic interaction present, the inhomogeneous term is generated rep- 
resented by T/(P) which is just real constant for a given discrete value of bound state mass 
P. So far all delta distributions in full p should be factorized and integrated out, thus we 
naturally take 

r ™ =r M + hJ*w$%hv (3 - 12) 

where p p (a, rf) is assumed to be a real function, and not a delta distribution. The nontrivial 
part of integral representation corresponds exactly with the function Y p , noting its structure 
is fully driven by pure triplet interaction of the Higgs. Furthermore, we explicitly choose 
n = 1 in the integral representation (13.121) . following the most easy integral representation 
of the inhomogeneous term in the expression, i.e. f k V p (k,p,P)G®(k,P). In fact, These 
integrals correspond to the scalar triangle Feynman diagram. Indeed, written more explicitly, 
we get 

/ V p (k,p,P)G [2] (k,P) = Ag 2 U2 I \D(k-p,M 2 ) + D(k + p,M 2 )]G [2] (k,P) 

Jk J k L J 

+ 36gj n ( \D(k-p,M 2 ) + D(k + p,M 2 )]G^(k,P) 

J k L -I 

= 8g 2 U2 I A112 (Pp) + 72g 2 in I A11 i(P,p) 

8#i 12 PAii2(a, z) + 72gf u p A111 (a, z 



-i 



(a - M?)F(ct,z;P,p) v 1 



where the functions p A s are listed in the Appendix A in Rel. (1A3j) . The factor (a — M 2 ) in 
the last line typically factorizes and here it is the matter of our convenient convention. 

From human effort point of view the choice n = 1 represents most convenient one, however 
recall here the known property of superrenormalized model studied yet: The function p^ 
is more smooth for larger n. The BSE were solved for generalized Wick-Cutkosky models 
for the lowest value n — 1, 2 in practice (the exception is the original Wick model, where n 
labels different energy states and is not arbitrary at all, this is because of massless boson). 
The later choice usually lead to meaningful improvement of numerical accuracy. In our case 
with n = 1 we expect 1-2 percentage accuracy in P 2 eigenvalue identification. 

In order to express the second term in (I3.5P we use the integral representation for T p , 
hence we can write 



1 oo 



J V p {k,p,P)G l2] (k,P)T p {k,P) = J dz Jdap 2 {a,z)x (3.14) 

— 1 — oo 

{4#ii2 [Itm(P,P\ °> z ) + 4iia(-P,P; a, z)] 

+ 3 %iii [Itm(P,p;a,z) + I uni (P,p;a, z)]} , 
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where 



I t n 2 (P,p;a,z) = f D(q + P/2, M 2 )D(-q + P/2, M 2 ) ffi q >*® , (3.15) 
Jk t [a, z; P, q) 

I u ii 2 (P,p;a,z) = [ D(q + P/2, M 2 )D(-q + P/2, M 2 ) ffi + ^2 
jfc -r (a, z; P, 

and the function lun, i«iii can be obtained from Itn 2 , luivz through the simple change 
M 2 -> Mi . 

The functions J t , 7 U satisfy the following integral representation 

1 oo . 

I t , u (P,P;a,z) = Jd V J da (3.16) 

where we have suppressed the propagators indices for brevity and where the functions p tiU 
are reviewed in the Appendix O As simple consequence of crossing symmetry we have 

p u (a,r],a,z) = p t (a,-T],a,z) . (3.17) 

At this point we can assume the validity of theorem of uniqueness derived for general 
PTIR which implies that the original BSE is equivalent to the solution of the following 
regular real integral equation: 

1 oo 



p p (a,r)) 



T I (P)p I {a, rj) + J dz J da p p (a, z)V(a,r],a,z) 



(3.18) 



a -Ml 
where 

Pi(<x,v) = 8g 2 112 pAii2(a,v) + 72gf u p A111 (a, 77) , 

V(a,r),a,z) = 4g 2 12 [p tll2 (a,r],a,z) + p tll2 (a, -r],a, z)] 

+ SQgf n [p tlll (a, 77, a, z) + p tlll (a, -77, a, z)] . (3. 19) 

As a consequence of the quartic constant interaction term Eq. (13. 18[) is inhomogeneous and 
necessarily coupled to the Eq. (13.91) . Using the integral representation for T p 

1 00 

Vj R J dz J dap p (a, z)Ip(P 2 ; a, z) 
rv(p) = =i — m , (3.20) 

where Vj R is the renormalized constant interaction and the suitable expression for the loop 
integrals If and Ib are reviewed in the Appendix [Bl 

In this work we use the momentum subtraction renormalization scheme with zero mo- 
mentum scale. Hence 2 ^ 

Vj - QX nn + p 2 -zr M - 2 + p^W 2 (3 ' 21) 

where Af xll is the renormalized quartic coupling of heavier Higgs mass eigenstates Hi, the 
others couplings stem from superrenormalizable interaction, hence they do not need to be 
renormalized at all. For a fixed bound state mass P 2 the function Tj(P) is a simple constant, 
affecting effectively the strength of the inhomogeneous term. Note that there is a peak for 
Pjf' (P 2 ) when approaching the continuum and the one loop effect can be quite large. 
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IV. RESULTS 



The inhomogeneous BSE requires renormalization since the function p p is linearly present 
at each term of BSE rhs. After a suitable normalization then the BSE for the weight function 
(I3.18p has been solved by the method of iterations. The coupling constants have been varied 
to reach the real discrete spectrum of Higgsonia. 

Firstly, let us mention the result for the SM Higgsonium as it follows from our method 
of BSE and compare with j^] . In this case, the only experimentally constrained input is 
the Hide's vev. v = 275/V2GeV while the Higgs mass and the cubic coupling depends on 
presently unknown coupling A, these satisfy = v^2At>, A3 = 2\v. As it is known, there are 
no bound states below a certain critical coupling of A3. The first state appears with the mass 
close to the on shell masses of Higgses M = 2m/j (M = 1.95m^ was taken in practice) and it 
is formed when rrih = 1.28TeV. This value is in a reasonable agreement with H[, however the 
input mass of such Higgs is ruled out by electroweak precision measurements. Furthermore, 
for such a fat Higgs, the Higgs sector of SM represents strongly coupled field theory and 
our BSE solution becomes only a rough estimate. In addition, being more realistic and 
switching on the top quark Yukawa, such fat Higgs becomes broad resonance and its fast 
decays should prevent the formation for bound states. The inhomogeneous BSE alone turns 
out to be not well controlled tool at all in this case. 

To achieve the numerical solution of Higgsonium in xSM we have firstly adjusted the 
couplings between the gauge eigenstate to get physical masses and couplings between the 
the mass eigenstates H\^Hi and then we look for the BSE solution varying the mass P 2 . 
A more convenient approach where the bound state mass has been fixed and some of the 
couplings have been varied has been followed as well. Both approaches have been found 
equivalent, however feasibility of getting the numerical solution can be different. The later 
approach, if needed, requires re-identification of the original gauge eigenstates couplings. 

As we can not scan the all parameter space, we choose some exemplatory group, we 
suppose is relevant for calculation of a bound states. An example of the parameters in 
used as an input is as the following: v = 275GeV, A = 0.20, 5\ = 1.20 * v, 62 = 0.40, 
K2 = 0.10 * v, K3 = 5.0 * v, K4 = 0.20. It leads to two massive eigenstates with Mi = 
179.5GeV, M 2 = 177.7GeV, and couplings g m ~ 2S0GeV,g 22 2 ^ 4Q0GeV , Ann ~ 0.33 
and the appropriate mixing is cosO = 0.696, such large value of mixing is physically justified 
since the masses M 1; M 2 are in the range allowed by electroweak precision test [3|. In this 
case we found the solution of BSE that produces the bound states 20 % lighter than the 
Higgs production threshold : 

M B = 1.6 x M x = 286GeV . (4.1) 

In general it is clear that the trilinear couplings must be large enough to ensure bound 
states, otherwise we get only continuum spectrum of free Higgses in a production experi- 
ments. As expected and observed, increasing trilinear couplings the mass of the bound state 
decrease and goes away from the sum of on-shell constituent masses. 

At given stage, being limited by numerical accuracy and convergence, we did not perform 
ultimate search for the limiting case P = 2Mui solution, neither we were looking for today 
physically not well motivated bootstrap -deeply bounded state solutions Mb « (Mi + 
M 2 ). Instead of, we simplify situation and we kept the masses of Higgses fixed and vary 
only trilinear couplings (7m, (7222 by a common prefactor C, #111, #222 — > Cgin, Cg222 m the 
parameter space of mass eigenvalue physical states. In this way we have found trilinear 
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couplings have changed by the significant prefactors C = 0.41 (gm = 114.8GeV, #222 = 
164GeV) in order to get the following bound state mass eigenvalue: 

M B = 0.95 * 2Mi = 95.53%(M! + M 2 ) , (4.2) 

which value is already rather closed to the threshold. In other words: we predict weakly 
bounded Higgsonia in the singlet extension of SM with natural trilinear couplings (measured 
in units of Higgs masses), while the deep bound states are possible only when the coupling are 
un-naturaly tuned to the artificial and perhaps un-natural large values. Increasing the mass 
Ms towards the threshold we still expect some substantial decrement of trilinear couplings, 
however the BSE did not provided stable and arbitrarily accurate solutions to the date. 

In general, the lighter Higgs we have and for the bound state mass closed to the thresholds 
M1+M2 we can expect that the presented form of the BSE is built up already in quite reliable 
approximations. Being not so far from the threshold, the next order irreducible diagrams 
(e.g. cross diagrams etc.) contribution is getting small since we stay effectively in the weak 
coupling regime (note that higher orders are suppressed by powers of l/(47r) 2 ). We do not 
present the BSE solutions for very deep bound state, which in principle can be interesting 
issue of pure numerical interest. For a very deeply binded Higgs we have no justification 
of approximations we made. Also, as we can learn from QCD-like studies, the importance 
of dressing constituent (via Schwinger- Dyson studies) could be addressed for the strong 
coupling systems. This issue is happily much less urgent for a weakly bounded Higgsonia, 
which as our study shows up, much likely exist in the singlet extension of SM. 



V. CONCLUSION 

Using a BSE formalism we performed a simple search for scalar bound states in the Higgs 
sector of SM and mainly in xSM extension of SM. In the second case, in order to have a 
reasonable model which is not completely ruled out by electroweak oblique correction con- 
straints we have assumed reasonably light scalars with large mixing. In this model, there 
are no light deeply bounded states unless the new cubic coupling is very large. On the other 
side, we predict bound state of the mass compared to the sum of the masses of the con- 
stituents, i.e. approximately non-relativistic bound states. This is the main conclusion from 
the numerical inspection of assumingly most promising regime of xSM parameter space. The 
Higgsonia productions can be expected in xSM with large mixing and in similar models with 
more singlets and doublets as well. Our fully relativistic calculation shows that Higgsonia 
can appear in simple realistic extensions of SM. If some sort of Higgsonia will appear realized 
in the nature, then the method presented here is improvable by considering improved kernels 
of the model. For instance, taking into account the selfenergy contributions or considering 
more complicated cross exchanges of Higgses is straightforward through the perturbative 



introduction of appropriately known integral representations of such kernels [111]. [16 
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APPENDIX A: INTEGRAL REPRESENATION FOR I A (P,p) 

In this appendix we review the integral representation for the following scalar triangle 
diagram 

I A (P,p)=ij -^D(q + P/2,M 2 )D(-q + P/2,Ml)D(q-p,M 2 ), (Al) 

This is just the scalar triangle diagram, hence the expression is finite and do not require 
any renormalization. The integral representation is easily derivable by following the Feyn- 
man parametrization tricks and we refer the reader to the Appendix G of the work [l5| for 
the details. 

The amplitude reads Ia 

1 oo 

W,P,*)=(*l* {a _$$l mPir y (A2) 



-1 — OO 



M= _J_ £ e W e(i- f ,)ep) 

32^ t± \a - *S - S\ 

i 

a + M$- Ml ± VP 
l± = 2fr=S) ' (A4) 

D = (a - Ml - Ml) 2 - AMl{Ml - S) , (A5) 

S = ^(l-z 2 ) (A6) 
Further note that, one can derive the following useful constraint 



where we remind 



a > a min = Ml + Ml + 2^M| (M 2 - S) . ( A7) 

This can be utilized when enforcing the convergence of the numerical solution. 



APPENDIX B: ONE LOOP BUBBLE I B AND I F 

The one bubble integral Ib reads 

I B (P 2 ) = i J** i D(q + P/2, Ml)D(-q + P/2, Ml) , (Bl) 

which is log divergent and so requires subtraction. In momentum subtraction scheme it 
defines the value of a renormalized quantity (coupling constant in our case) at the momentum 
scale p 2 = fi 2 . The renormalized one loop correction can be expressed as 

4W) = Hp 2 ) - W) 



dl-AMl/u P 2 -u 2 
^ m .. ^ , (B2) 



P 2 — UJ UJ — jl' 
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noting that the function Is is real for underthreshold momenta P 2 < AM 2 which is the 
relevant regime for our bound state study. 

The contact term leads to the following type of the integral in BSE 



I F (P 2 ; a,z)=ij ^ D(q + P/2, M 2 )D(-q + P/2, M 2 



(27T) 



F(a, z; P, q) 



(B3) 



This bubble-like contribution is finite since softened by the factor 1/F which stems from the 
bound state integral representation. The integral can be solved analytically by the manner 
described bellow. 

First let us use Feynman trick and match the propagator together, explicitly written we 
have 

If drj 



D (q + P/2, Ml)D { - q + P/2, Ml) = -J W(Ml m 9 _ p)]2 , 
matching further with the third fraction in the expression flB3j) we get 



D(q + P/2,M 2 )D(q-P/2,M 2 [ 
F(a, z; P, q) 



(B4) 



1 1 

dy / dt] 

J 1 [k.P[T] — z)y + (a — m 2 )y + k 2 + k.Pz + ^ — a] 3 



Integrating over the momenta we obtain 

1 

I F (P 2 -a,z) 



p2 

(B5) 



2(4vr) 2 (a 



b^w 2 1 dy I dr] i^z 

' -1 4 V 



j3 2 ) — m 2 y — a(l — y) + ie 



(B6) 



where (3 = (rj — z)y + z. After the substitution r] — > j3, the integration over the variable y 
and some trivial algebra we can get the following expression 



I F (P 2 ;a,z) 



1 



2(47r) 2 (a - Mf) 



J dp In ^ + J dplnj^- 

z l-z ' _1 1 + z 



U 



where we defined 



u 



M 2 -^(1-P 2 ) 
a -Mi 2 



(B7) 



:bs) 



APPENDIX C: INTEGRAL REPRESENTATION FOR I t AND I u 

The leading term in the BSE is a weighted sum over the following integral: 



I t (P,p; a, z) 



d 4 q 
(2^ 



D(q + P/2, M 2 )D(-q + P/2, M 2 ) 



D(p-q,M 2 ) 
F(a, z; P, q) 



(CI) 



The appropriate expression is derived in the paper [llj and [15|. Note, there are few misprints 
thorough the derivation in the original paper and for detailed derivation we recommend the 
later work 15 . 
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1 oo 



WW,*) = fdr, f da P t%;l\ a ' Z) P v (C2) 
J J (a - Mf)F(a,mP,v) 



-1 — oo 



a-M?)F(a,r);P,p) 



o(avaz)- 1 V > r e(x,(T))9(l-x,(T))e(D(T)) 

r), a, z) - -— — 2^2^ j ^ ' \^ 6 ) 

z l 47r ) T i=± \a - - S\ 

a + Ml-R{T)±VD 

X±(T) = ' (C4) 

D(T) = (a + Mj - R{T)f - AMj{a - S) , (C5) 



where 



R(T) = Ta + (1 — T)M\ , (C6) 
and the symbol J2t is shorthand notation for the following sum 

E /CO = /(o) - *(* - »j)/Cr + ) - % - *)/(r_) , (C7) 

where / is an arbitrary function and 

T ± = i±l (C8) 
± l±z V 7 

The term that is generated due to the u — term gives only a subleading contribution. Its 
derivation can be very easily proceed by considering the change t — > u. Explicitly we have 

I u (P,p;a,z)=i[ ^- A D( g + P/2, Mf)D(-q + P/2, M?) gjg + , ( C9 ) 

7 (2tt) 4 F(a,z;P,q) 

Hence, using the fact that I u (P,p; a, z) = It(P, —p; a, z) we get 



1 oo 



I u (P,p;a,z) = fdr, f da- ^( p v (CIO) 

J J (a - M{)F(a,r};P,-p) 

— 1 — oo 

which can be written as 



1 °° p (a —7] a z) 

— 1 — oo 

since P(a, 77; P,p) = F(a, —77; P, —p), or in other words p u (a, 77, a, 2;) = pt(a, —77, a, 2;). 
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